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Abstract 



o 

We studied the temperature and abundance distributions of intra-cluster medium (ICM) in the Abell 262 
■ cluster of galaxies observed with Suzaku. Abell 262 is a bright, nearby poor cluster with ICM temperature 

of ~ 2 keV, thus providing useful information about the connection of ICM properties between groups and 
clusters of galaxies. The observed spectrum of the central region was well-represented by two temperature 
models, and the spectra for the outer regions were described by single temperature model. With the XIS 
instrument, we directly measured not only Si, S, and Fe lines but also O and Mg lines and obtained those 
O ' abundances to an outer region of 0.43 rigp for the first time. We found steep gradients for Mg, Si, S, 

and Fe abundances, while O showed almost flat abundance distribution. Abundance ratios of a-elements 
^ ■ to Fe were found to be similar to those of the other clusters and groups. We calculated the number ratio 

of type II to type la supernovae for the ICM enrichment to be 3.0 ± 0.6 within O.lrigo, and the value was 
consistent with those for other clusters and groups. We also calculated metal mass-to-light ratios (MLRs) 
for Fe, O and Mg (IMLR, OMLR, MMLR) with B-band and K-band luminosities of the member galaxies 
^ , of Abell 262. The derived MLRs were comparable to those for other clusters with kT = 3-4 keV. 

Key words: galaxies: clusters: individual (Abell 262), galaxies: intergalactic medium, galaxies: 
abundances 



^ 1. Introduction 2001). Renzini (1997) and Makishima et al. (2001) sum- 
(yQ marized iron-mass-to-light ratios (IMLR) with B-band lu- 
. The metal abundances of Intra-cluster medium (ICM) minosity for various objects, as a function of their plasma 
^ ' carry a lot of information in understanding the chemi- temperature serving as a measure of the system richness, 
. ^ , cal history and evolution of clusters. X-ray observations and IMLRs in groups were found to be smaller than those 
' of ICM can immediately determine its temperature and in clusters. They also showed that the early-type galax- 
JH . metal abundances. Recent X-ray observations allow us to ies released large amount of metals which were probably 
. - - 1 measure temperature and metal abundance distributions formed through past supernovae explosions as shown ear- 
in the ICM based on the spatially resolved spectra. A Her by Arnaud et al. (1992). In order to obtain a correct 
large amount of metals of the ICM are mainly produced modeling of ICM, we need to know the correct temper- 
by supernovae (SNe) in early-type galaxies (Arnaud et al. ature and metal abundance profiles without biases (e.g., 
1992; Renzini et al. 1993), which are classified roughly as Buote 2000; Sanders & Fabian 2002). Especially for the 
type la (SNe la) and type II (SNe II). In order to know ICM of cooler systems, such as elliptical galaxies and 
how the ICM has been enriched, we need to measure the groups of galaxies, sensitive observations are required as 
amount and distribution of metals in the ICM. Because mentioned in Arimoto et al. (1997) and Matsushita et al. 
Si and Fe are both synthesized in SNe la and II, we need (2000). 

to know O and Mg abundances, which are synthesized Recent XMM-Newton and Chandra observations have 

predominantly in SNe II, in resolving the past metal en- enabled us to study the temperature structure, and de- 

richment process in ICM by supernovae. tailed properties of the heavy elements in the ICM were 

ASCA firstly showed the distribution of the heavy shown based on the high spatial resolutions and small 

elements, such as Si and Fe, in the ICM (Fukazawa point spread functions (Fabian et al. 2001; Fabian et al. 

et al. 1998; Fukazawa et al. 2000; Finoguenov et al. 2005; Fabian et al. 2006; Finoguenov et al. 2002; Xu et 

2000; Finoguenov et al. 2001), and BeppoSax also showed al. 2002; Gastaldello & Molendi 2002; Matsushita et al. 

the Fe distribution and its gradient (De Grandi & Molendi 2003; Matsushita et al. 2007b; Tamura et al. 2003; Sanders 
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& Fabian 2007; Humphrey & Buote 2006; Rasmussen & 
Penman 2007). In addition, the spatial distribution and 
elemental abundance pattern of the ICM metals were de- 
termined with the large effective area of XMM-Newton 
(Matsushita et al. 2007b; Tamura et al. 2004; Bohringer 
et al. 2005; O'SuUivan et al. 2005; Sanders & Fabian 
2006; de Plaa et al. 2006; de Plaa et al. 2007; Werner 
et al. 2006; Simionescu et al. 2008). However, the abun- 
dance measurements of O and Mg with XMM-Newton, 
in particular for the outer regions of clusters, are quite 
difficult due to the relatively high intrinsic background. 
Suzaku XIS can measure all the main elements from O to 
Fe, because it realizes lower background level and higher 
spectral sensitivity, especially below 1 keV (Koyama et 
al. 2007). Suzaku observations have shown the abun- 
dance profiles of O, Mg, Si, S, and Fe to the outer re- 
gions with good precision for several clusters (Matsushita 
et al. 2007a; Sato et al. 2007a; Sato et al. 2008a; Sato et 
al. 2008b; Tokoi et al. 2008). Comparing the Suzaku re- 
sults with supernova nucleosynthesis models, Sato et al. 
(2007b) showed the number ratios of SNe H to la to be 
3.5. 

Abell 262 (hereafter A 262) is a nearby cluster of galax- 
ies {z = 0.0163) characterized by a smooth and symmet- 
ric distribution of ICM, and is a richness class cluster. 
A 262 is located on a knot of the large-scale filament of the 
Pisces-Perseus supercluster. A 262 is a very useful system 
in studying the history of chemical evolution of groups 
and clusters, because the ICM temperature is about ~ 2 
keV which is intermediate between clusters and groups 
of galaxies. The central cD galaxy, NGC 708, is a host 
of a double-lobed radio source, B20149+35 (Parma et al. 
1986; Fanti et al. 1986). In the past, Einstein (David 
et al. 1993), ROSAT (David et al. 1996; Peres et al. 
1998; Neill et al. 2001), ASCA (White 2000; Fukazawa 
et al. 2004), XMM-Newton (Peterson et al. 2003), and 
Chandra (Blanton et al. 2004) observed A 262, and ASCA 
showed an average temperature of kT — 2.13 keV and 
an average metal abundance of 0.44 solar, respectively 
(Fukazawa et al. 2004). In a future paper (Gastaldello et 
al. 2008) we will investigate the mass and entropy radial 
profiles. 

This paper reports on results from Suzaku observa- 
tions of A 262 out to 27' ~ 540 h^Q kpc, correspond- 
ing to - 0.43 ri8o. We use Hq = 70 km s'^ Mpc~\ 
= 1 — = 0.73 in this paper. At a redshift of 
z — 0.0163, 1' correspon ds to 19.9 kp c, and the virial 
radius, rigo = 1.95 h^f^gy^k{T)/ 10 keV Mpc (Markevitch 
et al. 1998), is 1.25 Mpc (63') for an average tempera- 
ture of k{T) = 2.0 keV. Throughout this paper we adopt 
the Galactic hydrogen column density of A^h = 5.38 x 10^° 
cm~^ (Dickey & Lockman 1990) in the direction of A 262. 
Unless noted otherwise, the solar abundance table is given 
by Anders & Grevesse (1989), and the errors are in the 
90% confidence region for a single interesting parameter. 
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Fig. 1. Combined XIS image of central observation in the 
0.5—4.0 keV energy range. The observed XISO, 1, and 3 
images were added on the sky coordinate after removing 
each calibration source region, and smoothed with cr = 16 
pixel ~ 17" Gaussian. Estimated components of extragalac- 
tic X-ray background (CXB) and instrumental background 
(NXB) were subtracted, and the exposure was corrected, 
though vignetting was not corrected. The white circles show 
the extracted regions in the spectral fits. 

2. Observations and Data Reduction 

2.1. Observation 

Suzaku observed the central and offset regions of 
Abell 262 in August 2007 (PI: F. Gastaldello and K. 
Matsushita). The observation log is given in table 1, and 
the XIS image in 0.5-4 keV is shown in figure 1. We an- 
alyzed only the XIS data in this paper, although Suzaku 
observed the object with both XIS and HXD. The XIS 
instrument consists of three sets of X-ray CCDs (XIS 0, 

I, and 3). XIS 1 is a back-illuminated (BI) sensor, while 
XIS and 3 are front-illuminated (FI). The instrument 
was operated in the Normal clocking mode (8 s exposure 
per frame), with the standard 5 x 5 or 3 x 3 editing mode. 

2.2. Data Reduction 

We used version 2.0 processing data, and the analysis 
was performed with HEAsoft version 6.4.1 and XSPEC 

II. 3.2aj. Details of the analysis method are given in Sato 
et al. (2007a). Here we give just a brief description of the 
data reduction. The light curve of each sensor in the 0.3- 
10 keV range with a 16 s time bin was also examined in 
order to exclude periods with anomalous event rates which 
were greater or less than ±3cr around the mean to remove 
the charge exchange contamination Fujimoto ct al. (2007), 
while Suzaku data was little affected by the soft proton 
flare compared with XMM data. The exposure after the 
screening was essentially the same as that before screening 
in table 1, which indicated that the non X-ray background 
(NXB) was almost stable during the observation. Event 
screening with cut-off rigidity (COR) was not performed 
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Table 1. Suzaku Observation logs for Abell 262. 



Region 


Seq. No. 


Obs. date 


(RA, Dec)* 
J2000 


Exp. 

kscic 


After screening 
(BI/FI) kscc 


center 
offset 1 
offset2 


802001010 
802079010 
802080010 


2007-08-17X04:13:43 
2007-08-06X21:27:38 
2007-08-08X01:52:10 


(01'^52™46?1, -h36°09'33") 
(0l'^54™08?3, -h36°16'07") 
(01^52"14?4, -|-36°26'08") 


37.2 
54.6 
54.7 


37.0/37.0 
54.3/54.3 
54.3/54.4 



* Average pointing direction of the XIS, written in the RAJNFOM and DECJNTOM keywords of the event FITS files. 



in our data. In order to subtract the NXB and the extra- 
galactic cosmic X-ray background (CXB), we employed 

the dark Earth database by the "xisnxbgcn" Ftools task, 
and employed the CXB spectrum given by Kushino ct al. 
(2002). 

Wc generated two AnciUary Response Files (ARFs) for 
the spectrum of each annular sky region, A^' and , 
which respectively assumed uniform sky emission and 

1° X 1° size of the /3- model surface brightness profile, 
/? = 0.38 and = 0^47, in Fukazawa et al. (2004), by 
the "xissimarfgcn" ftools task (Ishisaki ct al. 2007). We 
also included the effect of the contaminations on the op- 
tical blocking filter of the XISs in the ARFs. Since the 
energy resolution also slowly degraded after the launch, 
due to radiation damage, this effect was included in the 
Redistribution Matrix File (RMF) by the "xisrmfgen" 
Ftools task. 

3. Temperature and Abundance Profiles 

3.1. Spectral Fit 

We extracted spectra from five annular regions of 0'-3', 
3'-6', and 6'-9' for the central observation, and 9'-17', 
17'-27' for the offset regions, centered on (RA, Dec) = 
(l'^52'"46?l, +36°09'33"). Xablc 2 lists the areas of the 
extracted regions (arcmin-^), fractional coverage of the an- 
nulus (%), the SOURCE_RATio_REG values (%; see caption 
for its definition) and the BI and FT coimts for the ob- 
served spectra and the estimated NXB and CXB spec- 
tra. Xhc fraction of the background, /bgd = (NXB -|- 
CXB)/OBS, was less than ^^50% even at the outermost 
annulus, although the Galactic component is not consid- 
ered here. Each annular spectrum is shown in figure 2. 
Xhe ionized Mg, Si, S, Fe lines are clearly seen in each 
region. Xhe O VII and O VIII lines were prominent in the 
outer rings, however, most of the O VII line was consid- 
ered to come from the local Galactic emission, and we 
dealt with those in the same way as Sato et al. (2007a) 
and Sato et al. (2008a). 

Xhe spectra with BI and FI for all regions were fitted 
simultaneously in the energy range 0.4-7.1 keV. In the 
simultaneous fit, the common Galactic emission compo- 
nent was included for all regions, while the A 262 emission 
was independently adjusted in each region. We excluded 
the narrow energy band around the Si K-edge (1.825- 
1.840 keV) because its response was not modeled correctly. 
Xhe energy range below 0.4 keV was also excluded because 
the C edge (0.284 keV) seen in the BI spectra could not be 
reproduced well in our data. Xhe range above 7.1 keV was 



also ignored because Ni line (~ 7.5 keV) in the background 
left a spurious feature after the NXB subtraction at large 
radii. In the simultaneous fits of BI and FI data, only 
the normalization parameter was allowed to take different 
values between them. 

It is important to estimate the Galactic component pre- 
cisely. Xhe Galactic component gives significant contri- 
bution especially in the outer regions, as shown in fig- 
ure 2. However, the ICM component is still dominant in 
almost all the energy range except for the O VII line. We 
assumed two temperature apec model (assiiming 1 solar 
abundance with redshift z = 0) for the Galactic compo- 
nent, and fitted the data with the following model for- 
mula: constant x [apec^^^i + apeCf^g^ + ICMcomponents), 
where we fixed the temperature of apec^^^i to be 0.074 
keV as given by in Lumb et al. (2002). As a result, the 
best-fit normalization of the cool (0.074 keV) component 
turned out to be 0. When we fix the two temperatures 
to be the value in Lumb et al. (2002), the resultant pa- 
rameters of the fit did not change significantly. Xhus, we 
concluded that the two temperature model of the Galactic 
emission (the cool component temperature was fixed to be 
0.074 keV) was enough to fit the A 262 data since the in- 
tensity of the Local Hot bubble in the A 262 direction 
seemed weak as shown in McCammon et al. (2002). In 
addition, we also examined the following model formula: 
constant x {apecf^g^+phabs x {apec^^^i + lCM component), 
however the values with this model were only a little 
worse than that with the above model. With this lat- 
ter model, the results for the ICM component did not 
change significantly. Xhe resultant normalizations of the 
apec models in table 3 are scaled so that they give the 
surface brightness in unit solid angle of arcmin^, and are 
constrained to give the same surface brightness and the 
same temperature for the simultaneous fits of all annuli. 

Xhe ICM spectrum for the central region, r < 3', was 
clearly better represented by two vapec models than one 
vapec model in the test. On the other hand, the ICM 
spectra for the outer regions, r > 3', were well-presented by 
single temperature model. Xhus, we carried out the simul- 
taneous fit with the following formula of the Galactic and 
ICM components: constant x {apec^g^i + apeci^g^+ phabs x 
(t;apeCo<r<27' + ''"^P^^o<r<3')- The fit results are shown in 
table 3. Xhe abundances were linked in the following way, 
Mg=Al, S=Ar=Ca, Fe=Ni, which gave good constraint 
especially for the offset regions. Xhe abundances were 
also linked between the two vapec components for r < 3' 
region. Results of the spectral fit for individual annuli 
are summarized in table 3 and figure 3, in which system- 
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Table 2. Area, coverage of whole annulus, SOURCE_RATIO_REG, and observed/estimated counts for each annular region. 



Kcgion 


Area ' 


Coverage 


^ SOURCE. 


* Energy 




BI counts ^ 






FI counts ^ 






( y/^Tn 1 n ^ 
1 cLiCllllli ) 










NXB 


CXB 


/BCD 




NXB 


CXB 


JBGD 


Center 


























0'-3' 


28.7 


100.0% 


13.6% 


0.4-7.1 


39,224 


316 


338 


1.7% 


57,652 


336 


534 


1.5% 


3-6 


84.8 


100.0% 


11.7% 


0.4-7.1 


37,674 


953 


871 


4.8% 


52,596 


970 


1,371 


4.5% 


6'-9' 


133.4 


94.3% 


9.8% 


0.4-7.1 


24,725 


1,551 


1,138 


10.9% 


33,040 


1,510 


1,705 


9.7% 


Offset 1 


























9'-17' 


121.2 


18.5% 


3.9% 


0.4-7.1 


13,956 


1,704 


1,448 


22.6% 


18,576 


1,922 


2,463 


23.6% 


17'-27' 


161.8 


11.7% 


3.3% 


0.4-7.1 


11,585 


2,718 


2,314 


43.4% 


13,507 


2,514 


3,398 


43.8% 


Offset2 


























9'-17' 


119.9 


18.3% 


4.6% 


0.4-7.1 


12,625 


1,945 


1,582 


27.9% 


15,062 


1,896 


2,268 


27.6% 


17'-27' 


162.6 


11.8% 


3.1% 


0.4-7.1 


9,884 


2,667 


2,244 


49.8% 


12,367 


2,640 


3,571 


50.2% 



SOURCE_RATIO_REG represents the flux ratio in the assumed spatial distribution on the sky (/3-model) inside the accumulation region to 
the entire model, and written in the header keyword of the calculated ARF response by "xissimarfgen" . 

* The inner three regions correspond to the central observation, and the outer two regions correspond to the offset observations, 
t The average values among four sensors are presented. 

t SOURCE.RATIO.REG = COVERAGE X J'J'°"' S{r) rdr/ S{r) rdr, where S{r) represents the assumed radial profile of Abell 262, and 
we defined S{r) in 1° X 1° region on the sky. 

^ OBS denotes the observed counts including NXB and CXB in 0.4—7.1 keV. NXB and CXB are the estimated counts. 
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Fig. 2. The panels show the observed spectra for the annular regions of A 262 which 
are denoted in the panels, and the data are plotted by red and black crosses for 
BI and FI, respectively. The estimated CXB and NXB components are subtracted, 
and the green and yellow lines show the best-fit model for the BI and FI spectra, 
respectively. The BI spectra of the ICM component are shown in blue and magenta 
lines, and the Galactic emission is indicated by cyan line. The energy range around 
the Si K-edge (1.825-1.840 keV) is ignored in the spectral fits. The lower panels show 
the fit residuals in units of a. 
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Table 3. Summary of the parameters of the fits to each annular spectrum of A 262. All annuli were simultaneously fitted. The 
values of Offset 1&2 correspond to the results of the summed spectra in the offset 1 and 2 regions. Errors are 90% confidence range 
of statistical errors, and do not include systematic errors. The solar abundance ratio of angr was assumed. These results are plotted 
in figure 3. 



Galactic 


Normi * 


(kcV) 


Norm2 * 


(keV) 








"•""-0.00 


0.074(fix) 


96+0-=''' 


n 1 74+0.036 

'^-0.043 






ICM 


Nornii 


kTi 
(keV) 


Norm2 ^ 


(keV) 


Normi / Norm2 


x7dof^ 


Center 














0'-3' 

3'-6' 
6'-9' 


365.4^1^7 
150.7111 
71.71}:? 


9 01+0.07 
^•"-•-^-0.06 
9 Q9+0.03 
^•"J^-O-Oa 
9 09+0.04 
^•'^^-0.04 


67.0l2i;0 


1 20+0-07 


K 41 + 1.71 

^•^-i-i.si 


895/497 

726/498 
648/498 


Offsetl 














9'~17' 
17' 27' 


c 9+0.2 
"•^-0.2 

1 7+0.1 

-■-•'-0.1 


9 -1 9+0.07 
^•-■^^-0.06 
1 fiQ+0.09 
J--Dy_0.04 








— 


364/382 
399/382 


Offset2 














9'-17' 

17'-27' 

total 


c 9+0.2 
^•^-0.2 


1 qfi+0.07 

^•'"-'-0.07 
1 4S+009 








— 


408/382 

530/379 
3970/3018 


Offsetl&2 














9'-27' 




1 94+0.04 
J--»^_0.04 








4149/3044 


ICM 





Ne 


Mg,Al 


Si 


S,Ar,Ca 


Fe,Ni 




(solar) 


(solar) 


(solar) 


(solar) 


(solar) 


(solar) 


Center 














0'-3' 
3'-6' 

6'~9' 


0.21^°;!° 

o.2otn^ 


o.88t°;f7 


0.771-J 
0.58iril 


0.84lH^ 
0.48t°;°^ 

o.52tn^ 


0.84^^? 
n fifi+o-07 

U.OO_Q 07 

n sn+o-09 

U.OU_o 09 


63+" "-^ 

^•O<^-0.03 

n 44+0-02 

U-*^-0.02 

36+0 03 
u.oo_o 03 


Offsetl 














9'-17' 
17'-27' 


42+0-20 

'-'•^^-0.20 
n 1 4+0-20 


A 77+0.28 
' ' -0.26 

n 42+0-27 

^•^^-0.14 


29+0-10 

U-^J^-0.17 


0.36;°;il 
n 1 's+o-io 


n 21+0-13 

^•^1-0. 14 

n 23+0-15 


-^0+0-04 

U.OU_o.o4 
-1 O+0.03 


Offset2 














9'-17' 

17'-27' 


n 24+0-^^ 


n S4+0-21 
"-'•''^-0.20 


o.o9l°:J^ 

0.29t0:l? 


n 23+0-10 

"•^^-0.10 

n 19+0-11 


(17+0-13 

U.U(_o.o7 

n 27+0-16 


17+0.03 
"-'•-l ' -0.03 

1 3+004 

^•-'^'J-0.03 


Offsetl&2 














9'-27' 


0.28triS 


0.54t°;ll 


0.24iro^ 


n 26+0-05 


n 1 7+0-07 

U.i ( _o.o7 


21+0-02 

"•■'-'--0.02 



* Normalization of the apec component divided by the solid angle, H", assumed in the uniform-sky ARF 
calculation (20' radius), Norm = J neUudV / {4tt{1 + z)^ D^) / CI" xlO-^" cm ^ axcmin ^, where D\ is 
the angular distance to the source, 
t The value is shown in Lumb ct al. (2002). 

t Normalization of the vapec component scaled with a factor of S0URCE_RATI0_REG / AREA in table 2, 
Norm = soL-R(-E_RATio_REG J n^^ndV / [4n{l+z)'^Dl] x lO-^o cm'^ arcmin'^, where Da is the angular 
distance to the source. 
§ All regions were fitted simultaneously. 
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Fig. 3. (a): Radial temperature profiles derived from the spectral fits for eacli annulus against the projected radius. Light-gray 
diamonds correspond to the results of the addition of offset 1 and 2 spectra. Black dashed lines show systematic change of the best-fit 
values of the center-offsetl region by varying the thickness of the OBF contaminant by ±10%. Light-gray dashed lines denote those 
when the estimated CXB and NXB levels are varied by ±10%. (b)-(f): Radial abundance profiles are plotted in the same way as in 
(a). 



atic error due to the OBF contamination and background 
(CXB -|- NXB) estimation are shown. 

3.2. Temperature Profile 

Radial temperature profile and the ratio of the vapec 
normalizations between the hot and cool ICM components 
are shown in figure 3(a) and table 3. The ICM tempera- 
ture of hot and cool components at the central region was 
~ 2.3 and ~ 1.2 keV, respectively, and the temperature 
decreased mildly to ^ 1.5 keV in the outermost region. 
Our results for the two temperature ICM model are con- 
sistent with the Chandra result (Blanton et al. 2004) for 
the central region within 3'. For the outer region (r > 3'), 
our results are also consistent with the previous Chandra 
result (VikhHnin et al. 2006), and the Chandra/XMM re- 
sults (Gastaldello et al. 2007). The radius of 27' - 540 kpc 
corresponds to ~ 0.43 rigo, and the temperature decline, 
observed in several other clusters, is clearly recognized in 
this system out to this radius. 

3.3. Abundance Profiles 

Metal abundances are determined for the six element 
groups individually as shown in figures 3(b)-(f). The four 
abundance values for Mg, Si, S, and Fe and their radial 
variation look quite similar to each other. The central 
abundances lie around ~ 1.0 solar, and they commonly 
decline to about 1 /5 of the central value in the outermost 
annulus. Note that, although the Fe abundance decreased 
to r >~ 20'; 0.3 rigo, the abundance at the outermost 



region was determined with Fe-L line as shown in figure 2 

because of the poor statistics around Fe-K line. On the 
other hand, the O profile looks flatter compared with the 
other elements. Because the results for the off'set regions 
had large errors, we examined the summed spectra for 
these regions. However, even if the spectra in the 9' < 
r < 17' and 17' < 27' regions for both the offsetl and 2 
observations were combined, the tendency of the radial 
abundance distributions did not change significantly. We 
also noted that, when we examined all abundances to be 
free in the fits, the resultant parameters did not change 
within the statistical errors. 

3.4- Systematic Errors and Uncertainties 

We performed the spectral fits with the vmekal model 
instead of the vapec model. As a result, the abundance 
profiles did not change within the statistical errors, al- 
though the Mg abundances with the vmekal model were 
~ 30-40% lower than those with the vapec model. 

We examined the systematic error of our results by 
changing the background normalization by ±10%, and the 
error range is plotted with light-gray dashed lines in fig- 
ure 3. The systematic error due to the background esti- 
mation is almost negligible. The other systematic error 
concerning the uncertainty in the OBF contaminant is 
shown by black dashed lines as shown in figure 3. A list 
of x^/dof by changing the systematic errors is presented 
in table 4. Although the fit is not statistically acceptable 
due mainly to the very high photon statistics compared 
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Table 4. List of x^/dof for the fits of tiie nominal and considering the systematic errors such as contaminant of OBF and background 
level. For details, see text. 



Region 


nominal 


contaminant 


back^ 


ground 






+ 10% -10% 


+10% 


-10% 


All 


. 3970/3018 


3769/3018 4350/3018 


4041/3018 


3997/3018 



with the systematic errors in the instrumental response, 
these results are useful to assess whether each element 
abundance is reasonably determined or not. We note that 
Ne abundance is not reliably determined due to an overlap 
with the strong and complex Fe-L line emissions, however 
we left these abundance to vary freely during the spectral 
fit. 

Oxygen abundance is strongly affected by the Galactic 
emission as shown in Sato et al. (2007a). However, we 
have shown that the use of different temperature models 
for the Galactic emission did not significantly affect the 
abundance results as mentioned in subsection 3.1. We also 
examined radial surface brightness profiles of the O VII 
and O VIII emission lines in order to see if they are spa- 
tially uniform or peaked at the cluster center. The surface 
brightness of the lines was derived by fitting the annular 
spectrum with a power-law + gaussian + gaussian model 
for an energy range between 500 and 700 eV. In this fit, 
we fixed the Gaussian cr to be 0, and allowed the energy 
center of the two Gaussians to vary within 555-585 eV 
or 630-660 eV for O VII or O VIII line, respectively. The 
derived line intensities are shown in figure 4. There is 
a clear excess in the OVIII intensity profile towards the 
cluster center, while O VII one is consistent to be con- 
stant as already noticed by Sato et al. (2007a); Sato et al. 
(2008a); Sato et al. (2008b). This is a clear evidence that 
the OVIII line is associated with the ICM itself, on the 
other hand, O VII line may mainly come from the Galaxy 
or from the interplanetary space. 

4. Discussion 

4-1. Metallicity Distribution in the ICM 

The present Suzaku observation of A 262 shows the 
abundance distribution of O, Mg, Si, S, and Fe out to a 
radius of 27' ~ 540 kpc, which corresponds to ^ 0.43 rigg, 
as shown in figure 3. The Ne abundance shows a large 
ambiguity due to the strong coupling with Fe-L lines. 
Distributions of Mg, Si, S, and Fe are similar to each 
other, while O profile in the outer region shows a large 
uncertainty. We plotted abundance ratios of O, Mg, Si, 
and S over Fe as a function of the projected radius in fig- 
ure 5. The ratios Mg/Fe, Si/Fe and S/Fe are consistent 
with a constant value around 1.5-2, while 0/Fe ratio for 
the innermost region (r < 2') is lower around 0.5. In ad- 
dition, the 0/Fe ratio suggests some increase with radius. 
The Fe abundance at the outermost region has a large 
uncertainty compared with that in the inner regions, be- 
cause the value in the outer region was determined by Fe-L 
line as shown in figure 2. We note that these abundance 
profiles are not deconvolved and averaged over the line of 
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Fig. 5. Comparison of the abundance ratio for A 262 (red) 
with results for the AWM 7 cluster (black diamonds; Sato 
ct al. 2008a) and NGC 507 group (blue crosses; Sato et al. 
2008b). Fe abundance and the O/Fe, Mg/Fe, Si/Fe, and S/Fe 
abundance ratios in solar units (Anders & Grevesse 1989) 
are plotted against the projected radius scaled by the virial 
radius, rigOi in a^U panels. 

sight. 

Recent Suzaku observations have presented abundance 
profiles in several other systems: an elliptical galaxy 
NGC 720 (Tawara et al. 2008) , a group of galaxies HCG 62 
(Tokoi ct al. 2008) and NGC 507 (Sato et al. 2008b), the 
Fornax cluster and NGC 1404 (Matsushita et al. 2007b), 
and a cluster of galaxies Abell 1060 (Sato et al. 2007a) 
and AWM 7 (Sato et al. 2008a). All systems show very 
similar value of Si/Fe ratio, to be 1-1.5. Mg/Fe ratio is 
slightly higher in HCG 62, Abell 1060 and AWM 7 than 
in other systems. We compare metal abundances of A 262 
with those of AWM 7 and NGC 507 as shown in figure 5. 
Because AWM 7 (- 3.5 keV) and NGC 507 (- 1.5 keV) 
host cD galaxies at the center and are similar to A 262 
in the X-ray morphology. A 262 may connect groups of 
galaxies and poor clusters regarding our understanding of 
the cluster evolution. Although the Fe abundance of A 262 
is slightly lower than that of the AWM 7 and NGC 507 
in the central region, the Fe profile to the outer region 
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Fig. 4. Intensities of (a) O VII and (b) O VIII lines at each annulus of A 262 in units of photons cm ^ s ^ sr ^. 



looks quite similar. Abundance ratios of 0/Fe, Mg/Fe, 
Si/Fe, and S/Fe are also quite similar between the three 
systems. Therefore, the abundance ratios show closer sim- 
ilarity than the absolute abundance values among different 
systems. The efficiency of the metal enrichment may de- 
pend on parameters such as age, star-formation efhciency, 
and contribution from cD galaxies. However, the relative 
contribution of SNe la and II and the process of metal 
mixing in the ICM seem to be quite similar among differ- 
ent clusters and groups. 

Tamura et al. (2004) reported abundance ratios for 
19 clusters studied with XMM-Newton, and the mean 
Si/Fe ratio in cool and medium temperature clusters with 
kT < 6 keV was ~ 1.4. This is consistent with the Suzaku 
results for groups and poor clusters including A 262. Their 
0/Fe ratio, ~ 0.6, in the cluster core also agrees with the 
Suzaku results including our A 262 case. Matsushita et al. 
(2003) and Matsushita et al. (2007a) reported abundance 
ratios for M87 and the Centaurus cluster, respectively, 
based on XMM-Newton observations. M87 showed Mg/0 
ratio to be ~ 1.3 in the central region, and the Centaurus 
cluster indicated 0/Fe and Si/Fe ratios within 8' to be 
~ 0.5 and ~ 1.0, respectively, and they were consistent 
with our results. Sanders & Fabian (2006) also showed the 
abundance ratios for the Centaurus cluster with Chandra 
and XMM-Newton, and the radial abundance ratios of 
0/Fe, Mg/Fe, and Si/Fe to be 0.5-1.0, - 1.0, - 1.2, 
respectively, were consistent with our results. 

4.2. Number Ratio of SNe II to SNe la 

In order to examine the SNe la and SNe II contribution 
to the ICM metals, the elemental abundance pattern of O, 
Mg, Si, S and Fe was examined for the region within 0.1 
and 0.27 rigo- The abundance ratios at r ~ 0.3 rigo and 
0.3 < r < 0.4 riso, with the latter corresponding to the out- 
ermost region, look similar within the statistic errors. The 
relative abundance ratios to Fc were fitted by a combina- 
tion of average SNe la and SNe 11 yields per supernova, 
as shown in figure 6. The fit parameters were chosen to 
be the integrated number of SNe la (iVja) and the number 
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Fig. 6. Fit results of number ratios of elements to Fe of 
A 262. Top and second panel show the abundance number 
ratios in solar unit within 0.27 rigo (black) and 0.1 rigo (red) 
regions, respectively. Blue dashed and solid lines in the second 
panel correspond to the contributions of SNe la (W7) and 
SNe II within 0.1 rigo, respectively. Ne (atomic number = 10) 
is excluded in the fit. Third panel shows the comparison with 
solar abundance table of Anders & Grevesse (1989) (light-gray 
line) and Lodders (2003) (dashed light-gray line). Forth panel 
indicates ratios of data points to the best-fit model. Bottom 
panel indicates fractions of the SNe la contribution to total 
metals in the best-fit model for each element, respectively. 
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Region SNc la Model TVia A^ii/A^ia xVdof 



<0.1ri8o 


W7 


4.4 ±0.4 X 10« 


3.0±0.6 


8.6/3 


<0.27ri8o 


W7 


2.7±0.9 X 10^ 


3.2 ±1.0 


4.1/3 


< 0.1 nso 


WDDl 


5.1 ± 0.4 X 10^ 


1.8 ±0.5 


27.2/3 


<0.27ri8o 


WDDl 


2.9 ±1.1 X 10^ 


2.1±0.9 


19.6/3 


< 0.1 nso 


WDD2 


4.2 ±0.4 X 10*^ 


2.9±0.6 


9.2/3 


<0.27ri8o 


WDD2 


2.5 ±0.9 X 10^ 


3.1 ±1.0 


7.6/3 



ratio of SNc II to SNe la (Nu/Ni^), because TVia could 
be well constrained by the relatively small errors in the 
Fe abundance. The SNe la and II yields were taken from 
Iwamoto et al. (1999) and Nomoto et al. (2006), respec- 
tively. We assumed a Salpeter IMF for stellar masses from 
10 to 50 Mq with the progenitor metallicity Z = 0.02 
for SNe II, and W7, WDDl or WDD2 models for SNe la. 
Table 5 and figure 6 summarizes the fit results. The num- 
ber ratios were better represented by the W7 SNe la yield 
model than by WDDl. The number ratio of SNe II to 
SNc la with W7 is ^ 3.2 within 0.3 riso, while the ratio 
with WDDl is ~ 2.1. The WDD2 model gave the resuh 
quite similar to those of W7. The resultant number ra- 
tios are consistent with those in Sato et al. (2007a). We 
also compared the abundance pattern of A 262 with the 
solar abundance. The third panel in figure 6 shows the 
comparison of the abundance pattern within r < 0.1 riso 
of A 262 with the solar abundance pattern in Anders & 
Grevesse (1989) and Lodders (2003). Our results of Mg, 
Si, and S fall between the abundance patterns of Anders 
& Grevesse (1989) and Lodders (2003). 

Almost - 80% of Fe and - 30% of Si and S were syn- 
thesized by SNe la in the W7 model, as demonstrated in 
the bottom panel of figure 6. The observed tendency of 

2 keV cluster is similar to those of ~ 3 keV clusters in 
Sato et al. (2007a). The values in table 5 imply that the 
A^ii/A^'ia ratio behaves in a similar manner for different su- 
pernova models between r < 0.1 riso and r > 0.1 riso- We 
note that the fit was not acceptable based on the value 
in table 5. As described in Sato et al. (2007a), the models 
adapted here (SNe yeild, Salpter IMF, etc.) are probably 
too simplified. 

4-3. Metal Mass-to-Light Ratio 

Combining the abmidance profile obtained with Suzaku 
and the 3-dimensional gas mass profile with Chandra 
and XMM-Newton (Gastaldello et al. 2007), we calcu- 
lated cumulative metal mass as shown in figure 7(a). In 
order to see whether the gas mass profile derived from 
Chandra/XMM observations could apply to the offset re- 
gions of Suzaku which were out of the field of view of 
Chandra/XMM observations, we compared three surface 
brightness profiles: 1. the Chandra/XMM observations 
which were used to derive the mass profiles, 2. the ASCA 
observations which were used to calculate ARFs, and 
3. the present Suzaku observation. However, we could 
not compare these three profiles directly because of the 
poor angular resolution and vignetting effect of Suzaku 



as shown in Sato et al. (2007a). Thus, we need to sim- 
ulate the surface brightness profiles using the parameters 
of the XMM/Chandra and ASCA observations with "xis- 
sim" Ftools task. As a result, those profiles were consis- 
tent with that of this Suzaku observation within r < 13', or 
260 kpc. On the other hand, the profile in r > 20' looked 
slightly diff'erent, and the observed Suzaku profile was 
slightly steeper than those of Chandra/XMM and ASCA. 
Thus, the derived metal mass has a systematic error by 
about a factor of 2. The derived iron, oxygen, and magne- 
sium mass within the 3-dimensional radius of r < 540 kpc 
are 3.1 x 10^, 2.4 x 10^°, and 1.6 x 10^ Mq, respectively. 
Errors of the metal mass plotted in figure 7(a) are taken 
from the statistical errors of each elemental abundance in 
the spectral fits, because the statistical errors of the fits 
are lager than those of the gas mass profiles in Gastaldello 
et al. (2007). 

We examined metal mass-to-light ratios for oxygen, 
iron, and magnesium (OMLR, IMLR, and MMLR, re- 
spectively) to compare the ICM metal distribution with 
the stellar mass profile. Although, historically, B-band 
luminosity has been used for the estimation of the stel- 
lar mass (Makishima et al. 2001), we calculated the K- 
band luminosity in A 262 based on the Two Micron All 
Sky Survey ^ with the Galactic extinction, Ak = 0.032, 
from NASA/IPAC Extragalactic Database (NED) in the 
direction of A 262. We used 2° x 2° data set centered at 
the A 262 position as shown in table 1, and subtracted 
the r > 1° region as the background. In addition, we de- 
projected the luminosity profile as a function of radius. 
The resultant luminosity within this Suzaku observation, 
r < 27', is 6.7 x 10^^ Lq in K-band. The radial luminosity 
profile is not directly plotted, but it can be inferred from 
the ratio of the metal mass in figure 7(a) divided by the 
metal mass to light ratio in (b). We calculated the inte- 
grated values of OMLR, IMLR, and MMLR with K-band 
within r ^ 540 kpc as shown in figure 7(b) and table 6, and 
their values turned out to be ~ 3.6 x 10~^, ~ 4.6 x 10~^, 
and ~ 2.4 x 10~^ Mq/Lq, respectively. The errors are 
only based on the statistical errors of metal abundance 
in the spectral fit, and the uncertainties of the gas mass 
profile and the luminosity of member galaxies are not in- 
cluded. Note that the derived luminosity profiles were 
not corrected the Suzaku PSF effect, because the errors 
of MLRs were dominated by the error of the metal mass 
compared with that of the light luminosity. 

^ The database address: http://www.ipac.caltech.edu/2mass/ 
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Integrated mass-to- light ratios of O, Mg, and Fe (OMLR, MMLR, IMLR ) with K-band luminosity in units of Mq/L^. 



Region (kpc/nso) OMLR MMLR IMLR 



<59.8/0.05 


4 5+^-^ 


X 10" 


4 


1 1+"-^ 


X 10" 


-4 


2 6+'^-^ 

^•"-0.1 


X 10" 


-4 


<119.5/0.10 


^•^-0.7 


X 10- 


3 


o 4+0.5 


X 10- 


-4 


7 4+0.3 


X 10- 


-4 


<179. 3/0.14 


4.2tn 


X 10" 


3 


7 4+1-0 

'•^-1.0 


X 10" 


-4 




X 10" 


-3 


<338.6/0.27 




X 10" 


2 


1 9+0.6 
■•-•^-0.5 


X 10" 


-3 




X 10" 


-3 


<537.8/0.43 




X 10" 


2 


^•^-1.3 


X 10" 


-3 


4.61°;? 


X 10" 


-3 
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Fig. 7. (a) Cumulative mass, M (< R), within the 3-dimensional radius, R, for O, Fc, and Mg in A 262, based on the combination of 
the abundance determination with Suzaku and the gas mass profile with Chandra and XMM-Ncwton. (b) Ratio of the O, Fe, and Mg 
mass in units of Mq to the K band luminosity in units of Lq (OMLR, IMLR, and MMLR, respectively) against the 3-dimensional 
radius. The black crosses, black diamonds, and gray diamonds show the OMLR, IMLR, MMLR, respectively. 



Table 7. Comparison of IMLR, OMLR and MMLR with B-band luminosity for all systems. 





IMLR 


OMLR 


MMLR 


r (kpc/ri8o) 


k{T) 


Reference 


Suzaku 














NGC 720 .... 


1 X lO-'' 


4 X 10-4 




25/0.04 


- 0.56 keV 


Tawara et aL (2008) 


NGC 5044 . . . 


2.6 X 10-3 


6.6 X 10-3 


1.6 X 10-3 


88/0.10 


~ 1.0 keV 






3.6 X 10-3 


9.4 X 10-3 


2.6 X 10-3 


260/0.30 




Komiyama et aL (2008) 




6.0 X 10-4 






327/0.38 




Buote ct aL (2004) 


Fornax 


4 X 10"'' 


2 X 10-3 




130/0.13 


~ 1.3 keV 


Matsushita et al. (2007b) 


NGC 507 ... . 


6.0 X 10-4 


2.6 X 10-3 


3.7 X 10-4 


120/0.11 


~ 1.5 keV 






1.7 X 10-3 


6.6 X 10-3 


1.1 X 10-3 


260/0.24 




Sato et aL (2008b) 


HCG 62 


2.0 X 10-3 


6.4 X 10-3 


1.0 X 10-3 


120/0.11 


- 1.5 keV 






4.6 X 10-3 


3.8 X 10-2 


1.5 X 10-3 


230/0.21 




Tokoi et aL (2008) 


A 262 


3.6 X 10-3 


1.2 X 10-2 


1.6 X 10-3 


130/0.10 


~2 keV 






6.7 X 10-3 


3.7 X 10-2 


2.7 X 10-3 


340/0.27 




This work 


A 1060 


5.7 X 10-3 


4.3 X 10-2 


2.4 X 10-3 


180/0.12 


~3 keV 






4.0 X 10-3 


4.3 X 10-2 


1.6 X 10-3 


380/0.25 




Sato et al. (2007a) 


AWM 7 


4.8 X 10-3 


2.6 X 10-2 


3.4 X 10-3 


180/0.11 


~ 3.5 keV 






7.6 X 10-3 


3.1 X 10-2 


6.7 X 10-3 


360/0.22 




Sato et al. (2008a) 


XMM-N(nvton 
CcuUiunis . . . 


4 X 10-' 


3 X U)-- 




190/0.11 


- 4 kvY 


ALusushil a ct ul, (2007a) 
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Fig. 8. Comparison of IMLR (a) and (c), and OMLR (b) and (d) with B-band luminosity to the other clusters and groups within 
~0.1 ri8o and ~0.25 riso region, respectively. 



Wc also derived the MLRs with B-band himinosity us- 
ing an appropriate color B — K = 4.2 for early-typo galax- 
ies given by Lin & Mohr (2004), along with the Galactic 
extinction, Ab = 0.267, from NED in the direction of 
A 262. The integrated values of OMLR, IMLR, and 
MMLR with B-band within r ^ 540 kpc turned out to 
be ~ 7.1 X 10-2, ~ 9.0 x lO'^, and ~ 4.7 x lO^^ Mq/Lq, 
respectively. Comparing these B-band MLRs with those 
of other clusters and groups measured within ~ 0.1 rigo 
and ~ 0.25 rigo, the A 262 points fall between those of 
clusters and groups, in particular within ~ 0.1 riso- In 
other radii, the A262 results are consistent with those for 
the poor clusters as shown in table 7 and figure 8. The 
tendency that smaller systems with lower gas tempera- 
ture tend to exhibit smaller MLRs, as shown previously 
by Makishima et al. (2001) and by Sato et al. (2008b), 
is clearly indicated by the present addition of the A 262 
data. 

5. Summary and conclusion 

Based on the Suzaku observation of A 262, we studied 
spatial distribution of temperature and metal abundances 
for O, Mg, Si, S, and Fe up to ~ 0.43 nso- The ICM 
temperature decreases mildly from 2.3 keV to ~ 1.5 keV 



in the outer region. The abundances of Mg, Si, S, and 
Fe drop from subsolar levels at the center to ~ 1/5 solar 
in the outermost region, while the O abundance shows 
a flatter distribution around ~ 0.5 solar compared with 
the other metals. The abundance ratios, 0/Fe, Mg/Fe, 
Si/Fe, and S/Fe for A 262 are similar to those of other 
clusters with kT = 3-4 keV and groups with kT ~ 1.5 
keV. The number ratio of SNe II to la which contributed 
for the ICM enrichment is 3.0 ±0.6, consistent with the 
values for other clusters and groups. The derived MLRs 
with B-band luminosity fall between the data for clusters 
and groups measured within a radius of ~ 0.1 rigo- Thus, 
the A 262 data gives us a good connection concerning the 
global picture of the metal enrichment process between 
clusters and groups of galaxies. 
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